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Abstract 

t" — We would like to show the existence of finite energy SU(2) Yang-Mills- 

Higgs particles of one-half topological charge. The magnetic fields of these 

^3 solutions at spatial infinity correspond to the magnetic field of a positive 

one-half magnetic monopole located at the origin and a semi-infinite Dirac 
string which carries a magnetic flux of ^ going into the center of the sphere 
r~| at infinity. Hence the net magnetic charge of the configuration is zero. The 

solutions possess gauge potentials that are singular along one-half of the 
z-axis, elsewhere they are regular. There are two distinct configurations of 
these particles with different total energies and magnetic dipole moments. 
Their total energies are found to increase with the strength of the Higgs 
field self-coupling constant A. 
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1 Introduction 

The SU(2) Yang-Mills-Higgs (YMH) field theory in 3 + 1 dimensions, with the 
Higgs field in the adjoint representation possess magnetic monopole configurations 
[T]-jl]. These monopole solution belongs to the category of solutions which are 
invariant under a U(l) subgroup of the local SU(2) gauge group. The 't Hooft- 
Polyakov monopole solution with non-zero Higgs mass and self-interaction belongs 
to such category and it is the first finite energy monopole solution. This numerical 
monopole solution of unit magnetic charge is spherically symmetric pQ. How- 
ever, finite energy monopole configurations with magnetic charges greater than 
unity pE] cannot possess spherical symmetry [5]. Exact monopole solutions ex- 
ist only in the Bogomol'nyi-Prasad-Sommerfield (BPS) limit [3]-[I]. Outside this 
limit and when the Higgs field potential is non-vanishing, only numerical solu- 
tions are known. Other exact solutions found are the A-M-A (antimonopole- 
monopole-antimonopole) and vortex ring solutions and various other mirror sym- 
metric monopole configurations |6j. Numerical BPS monopole solutions with no 
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rotational symmetry have also been discussed [7j. Other numerical finite en- 
ergy monopole solutions include the MAP (monopole-antimonopole-pair), MAC 
(monopole-antimonopole-chain), and vortex ring solutions [8]. 

However, most of these monopole solutions reported so far are of integer topo- 
logical monopole charges. Hence it is our purpose in this paper to discuss the 
existence of finite energy one-half monopole magnetic charge solutions. These 
one-half monopole solutions are numerical solutions that are solved by first ob- 
taining the exact one-half monopole solutions at large distances and then fixing 
the boundary conditions at small distances and along the z-axis. 

We analyse four types of one-half monopole solutions which we label as Type 
Al, A2, Bl, and B2. Our analysis shows that the Type B (Bl and B2) solutions 
are exactly 180° rotation of the z-axis about the origin, r = 0, of the Type A 
(Al and A2) solutions (respectively). However for a particular value of Higgs self- 
coupling constant, A, the Type 1 (Al and Bl) and Type 2 (A2 and B2) solutions 
possess different total energies and magnetic dipole moments. Comparisons of the 
gauge invariant form of the gauge field potentials and Higgs fields are also made 
between the four one-half monopole solutions to confirm our conclusion that the 
Type 1 and Type 2 solutions are different one-half monopoles. 

The gauge potentials of the Type A (Type B) solutions are singular along the 
negative (positive) z-axis. The 't Hooft magnetic fields of these solutions at spatial 
infinity correspond to the magnetic field of a positive one-half monopole located 
at the origin, r = 0, and a semi- infinite Dirac string located on one-half of the 
z-axis which carries a magnetic flux of y going into the center of the sphere at 
infinity. Hence the net magnetic charge of the configuration is zero. The location 
of the Dirac string is always on the opposite side of the z-axis of its energy density 
concentration. 

When A = 0, the dimensionless total energy of the Type 1 one-half monopole 
solutions is 0.509 and the that of the Type 2 solutions is 0.525. Hence the total 
energies of these one-half monopole solutions seem to be higher than the dimen- 
sionless total energy of the BPS solution which is exactly |. The total energies of 
these one-half monopole solutions are calculated for various strength of the Higgs 
field self-coupling constant and they are found to increase with A. 

Recently we have constructed one-half monopole configurations which are ax- 
ially symmetric and possess a Dirac string along one-half of the z-axis [9]. The 
magnetic fields are the usual magnetic field of a one-half monopole charge that is 
spherically symmetric and radial in direction. 

Other works on one-half monopoles in the Yang-Mills (YM) theory include the 
work of Harikumar et al. pH]. They demonstrated the existence of generic smooth 
YM potentials of one-half monopoles. However no exact or numerical solutions 
have been given. Exact one-half monopole axially symmetric solutions and one- 
half monopole mirror symmetric solutions with Dirac-like string have also been 
discussed in Ref. JTT] . 

We briefly review the SU(2) YMH field theory in the next section. The mag- 
netic ansatz used in obtaining the one-half monopole solutions and some of its 
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basic properties are discussed in section 3. We present the method used in ob- 
taining the exact one-half monopole solutions at large distances in section 4 [9]. 
The numerical finite energy one-half magnetic monopole solutions and some of its 
properties are presented and discussed in section 5. We end with some comments 
in section 6. 

2 The SU(2) YMH Theory 

The SU(2) YMH Lagrangian in 3+1 dimensions with non vanishing Higgs potential 
is given by 

11 1 2 

C = pa paiiv _ _ D H$* D _ _A($ a $ a - ^-f . (1) 

4 2 " 4 v A 

Here the Higgs field mass is /i and the strength of the Higgs potential is A which 
are constants. The vacuum expectation value of the Higgs field is £ = fi/y/X. 
The Lagrangian ([I]) is gauge invariant under the set of independent local SU(2) 
transformations at each space-time point. The covariant derivative of the Higgs 
field and the gauge field strength tensor are given respectively by 

D^ a = d^ a + ge abc A b ^ c , 
F; v = d^-d^ + ge^AlAl, (2) 

where g is the gauge field coupling constant. The metric used is g^ v = ( — h ++). 
The SU(2) internal group indices a,b,c = 1,2,3 and the space-time indices are 
fi, z/, a = 0, 1, 2, and 3 in Minkowski space. The equations of motion that follow 
from the Lagrangian ([T| are 

D^F^ = d"F^ + ge abc A b »F c ^ u = ge abc <$> b D„<l> c , 
D»D^ a = A$ a ($ b $ fe -£ 2 ). (3) 

In the limit of vanishing fi and A, the Higgs potential vanishes and self-dual so- 
lutions can be obtained by solving the first order partial differential Bogomol'nyi 
equation, 

Bf ± A$ a = 0, where Bf = -\e ijk F« k . (4) 

The electromagnetic field tensor proposed by 't Hooft pQ upon symmetry break- 
ing is 

F^y = $> a F^--e abc &D fl $ b D u i> c , 

= d fl A u -d u A tM --e abc ^ a d tM $%$ c , where (5) 
9 

= d^A v -d v A^ andH tll/ = --e abc &d^ b d u &, (6) 
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are the gauge part and the Higgs part of the electromagnetic field respectively. Here 
= ^ a A^, the Higgs unit vector, $ a = $ a /|$|, and the Higgs field magnitude 
|$| = y/ Hence the decomposed magnetic field is 

Bi = --tijkF jk = Bf + Bf, (7) 

where Bf and Bf are the gauge part and Higgs part of the magnetic field respec- 
tively. The net magnetic charge of the system is 

M = — [ d { Bi d 3 x = — ld 2 Oi B^ (8) 



4tt J 4vr 
The topological magnetic current 



K = ±- e, upa e abc d^ a d^ b d°&, 
Sir 

is also the topological current density of the system. Hence the corresponding 
conserved topological magnetic charge is 

M H = 1 -Jd 3 xk = ^-J e l]k e ahc d % (^> a d^ b d k ^ d 3 x 



d 2 a l -e m e abc $ a d^ b d k V 



1 /* ,o / 1 

«7T J \g 

= Ifd^Bf. (9) 

The magnetic charge Mh is the total magnetic charge of the system if and only if 
the gauge field is non singular [13]. If the gauge field is singular and carries Dirac 
string monopoles Mq, then the total magnetic charge of the system is 

M = M G + M H , where 
M G = j d 2 a i e ijk {d j A k - d k Aj) 

= lJd 2 a lB f. (10) 

In the electrically neutral BPS limit when the Higgs potential vanishes, the 
energy is a minimum, 



E mm = T J d l (B a i ^ a ) d 3 x + J ^(Bf ± D^ a ) 2 d 3 x 

= T J d t {B^ a ) d 3 x = ——Mjjj (11) 

when the vacuum expectation value of the Higgs field is non zero. Hence the 
dimensionless minimum total energy is M#. 
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For a non BPS solution, its energy must be greater than that given by Eq. 
(11). The dimensionless value is given by 

E = — [ {B?B? + A$ a A$ a + -($ a $ a - £ 2 ) 2 } d 3 x > M H . (12) 
871-4 J 2 

The gauge potential A® and the Higgs field $ a can be gauge transformed by 
local SU(2) gauge transformations which can be written in the 2x2 matrix form, 

lo(x) = exp (^-a a h a (x)f(x)*j = cos Q/(^)J + iv a n a {x) sin Q/(^)J , (13) 

where n a (x) is a unit vector and a a is the Pauli matrix. The transformed gauge 
potential and Higgs field then takes the form, 

A 1 ; = cos fA«+ sin fe abc Alh c + 2 sin 2 U«(n b A b J 

+ {n a dj + sin fd^n a + 2 sin 2 f -t abc {d ^h b )n c ^ (14) 
$' a = cos /$ a + sin fe abc $ b h c + 2 sin 2 {h a (n b a b ). (15) 



3 The Magnetic Ansatz 

The magnetic ansatz j8] is given by 

gA* = -^M)n^ + — ^AM)n^ 

r v rsmu r v rsmu 

gA a = 0, g& = $ 1 {r,0)hi + $ 2 (r,e)fil (16) 

where Pi(r,9) = sin 9 ^2 (?",#) and P2(r,9) = sin 9 R,2(r,9). The spatial spherical 
coordinate orthonormal unit vectors are 

fj = sin^ cos0 5n + sin# sin0 5 i2 + cos# 8^, 
§i = cos 9 cos <p 5n + cos 9 sin0 8i2 — sin9 8^ 
4>i = - sin0 8 a + cos0 5 i2 , (17) 

and the isospin coordinate orthonormal unit vectors are 

h a r = sin9 cos ncj) 81 + sin 9 sinn0 8% + cos6> 8%, 
fig = cos 9 cos ncj) 81 + cos 9 sinn0 8% — sin# 8%, 
fA = — sin n<fr 8^ + cos n(j) 8%] where n > 1. (18) 

The 0-winding number n is a natural number. In our work here on one-half 



monopole we take n = 1. The magnetic ansatz (16) is form invariant under the 
gauge transformation 

u = exp (^cr a n^f(r, 9)*j , a a = Pauli matrices (19) 



and the transformed gauge potential and Higgs field take the form, 
1 



gA'i 



-{^ - defyhpi 

{Pi cos / + P 2 sin / + n[sin 9 - sin(/ + 9)]} h a 9 c 



1 



+ 

r sm( 

+ -{Rx + rd r f}h%fi 



1 



gA 

g<S> 



r sin 9 
,a — o 

o — u ' 
la 



{P 2 cos / - Pi sin / - n[cos9 - cos(/ + 6>)]}n"<; 



($i cos / + $2 sin /) n a r + ($ 2 cos / - $i sin /)n^. 



(20) 



The general Higgs fields in the spherical and the rectangular coordinate systems 



are 



g $ a = Q^x) fi a r + $ 2 {x)n a e + Q*{x)n% 

= 5 al + $ 2 (x) 5 a2 + $ 3 (x) 5 a3 , 



(21) 



respectively, where 



$i = sin^cosn0 $i + cos^cosn0 $2 ~ sinra^ $3 = |$| sinacos/3 
$ 2 = sin 9 sin ncj) $i + cos # sin n0 $2 + cos n0 $3 = | $ | sin a sin /3 
i> 3 = cos# $1 - sin6> $ 2 = |$| cosa. (22) 

The axially symmetric Higgs unit vector in the rectangular coordinate system is 

$ a = sin a cos /3 5 al + sin a sin /3 5 a2 + cosa 5 a3 , (23) 



cos« = g(r, 9) cos 9 — h(r, 9) sin 9, (3 = nc, 
9(r,e) = t^t, fc(r,0) - 



(24) 



m 



From Eq. (23) and (24), the Higgs field (16) and the gauge transformed Higgs 



field (20) can be written in terms of a as 



|$(r,0)|(cos(a - 6) h a r + sin(a - 9)h a 9 ), 
|$(r,0)|(cos(a' - 9) h a r + sin(a' - 0)ng) 



(25) 



where a' = a — f. 

Similarly the Higgs part of the 't Hooft mag netic field Q can be reduced to 

g Bf = -ne ijk d j cos ad k (f). (26) 

The gauge part of the mag netic field Q can be written in similar form 

gBf = -ne ijk dj cos « <9 fc 0, (27) 
1 



where, cos k 



7? 



(h(r,9)P 1 -g(r,9)P 2 



Hence the 't Hooft's magnetic field which is the sum of the Higgs part (26) and 
the gauge part (27) is given by 



gBi = -ne ijk d j (cosa + cosk) d k (p = -e ijk djA k , (28) 

where Ai is the 't Hooft's gauge potential. The magnetic field lines of the configu- 
ration can be plotted by drawing the contour lines of (cos a + cos k) = constant on 
the vertical plane = 0. The orientation of the magnetic field can also be plotted 
by using the vector plot of the magnetic field unit vector, 

- —de(cos a + cos K)fi + rd r (cos a + cos K)9i 

r>i = — (29) 
A/[r<9 r (cosa + cos ft)] 2 + [de(cosa + cos ft)] 2 

At spatial infinity in the Higgs vacuum, all the non-Abelian components of the 
gauge potential vanish and the non-Abelian electromagnetic field tends to 

CLoc = {d»A u -d u A,-- g e cd ^d,¥d u &W 

= F^ a , (30) 

where F^ is just the 't Hooft electromagnetic field. However there is no unique 
way of representing the Abelian electromagnetic field in the region of the monopole 
outside the Higgs vacuum at finite values of r [15J. One proposal was given by 't 
Hooft as in Eq. ([5]). The 't Hooft's magnetic field has the special property, that 
the magnetic charge density vanishes, that is d l B i = 0, when |$| ^ 0. However 
when |$| = 0, d l Bi ^ and the magnetic charges are located at these points. 
Hence with 't Hooft's definition of the electromagnetic field, the magnetic charges 
are discrete and reside at the point zeros of the Higgs field. Since this definition 
gives a discrete magnetic charge at a particular point, there is no magnetic charge 
distribution over space. 

Another proposal which is less singular was given by Bogomol'nyi j3] and Fad- 
deev [16] . 

Bi = B *(l)' S ' = Et (j)' (31) 



With this definition of the electromagnetic field (31), there will be a magnetic 
charge density distribution contributed by the non-Abelian components of the 
gauge field in the finite r region. A 3D surface plot of the magnetic charge den- 
sity distribution will be useful in determining the sign of the magnetic charge 
of the solutions. In the Higgs vacuum at spatial infinity, both definitions of the 
electromagnetic field ^ and (31) become similar. 
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4 The Exact Asymptotic Solutions 

The gauge transformation ( Jl9| , when applied on to the solution, 

gA a , gA« = A(r,9)5 a 3 &, g$ a = ^, (32) 

i a i n \ a cos 9 + b— I . 

where A(r, 0) = , a,o = constants. 

r sin 9 



will rotate the isospin direction, 5%, into the magnetic anstaz (16) and the trans- 
formed gauge potentials and Higgs field then take the form, 

gAT = - l -{-d e f^e)}hie i+ l -{d T f{T,9)}hih 

ft 

+ -^{sin0-sin(/(r,0) +6)[1 + rA(r,6) sin0] }h a e <j) i 

ft 

{cos 9- cos(/(r, 9) + 9)[1 + rA(r, 9) sin0]} 



r sin^ 



Ja 



g& a = £{™s(f(r,9)+9)h a r -sm(f(r,9) + 9)n a e }. (33) 

We also note that the angle a = —f(r, 9). Hence 

cos a = cos f(r, 9), cos k = — cos f(r, 9) + {1 + r sin 9 Air, 9)} , (34) 

and the 't Hooft gauge potential and net magnetic field of the monopole configu- 
ration are given by 

a cos 9 + b 



A.% 



r sin 9 



Tid 

gBi = -e ijk d j {l + rsm9A(r,9)}d k (l)= —n, (35) 
respectively. Hence one-half magnetic monopole solutions of the magnetic ansatz 



(16) can be obtained when n = 1 and a = 1/2. 

The Type Al solution is obtained by letting fir, 9) = —\9, a = |, b = |, and 
n = 1. The profile functions of the gauge potentials are given by, 

ipi = -, Pi = sm9 - ^sin^6»(l + cos6»), 

R x = o, P 2 = cos 9 - i cos ^0(1 + cos 9), 

3>i = £cos^, $ 2 = _£ sin Ifl. ( 36 ) 

The Type Al solution is obtained by letting f(r, 9) = \9, a = |, b = |, and n — 1. 
The profile functions of the gauge potentials are given by, 

1 13 
i>i = -Pi = sin0 - - sin -0(1 + cos6»), 

1 3 

R x = 0, P 2 = cos 9 - - cos -9(1 + cos0), 

$i = fcos^fl, $ 2 = -£sin^0. (37) 



Both Type A solutions possess a magnetic charge of ^ at r = and 't Hooft's 

gauge potential, Ai = { 2°^ n e } wn i cn is singular along the positive z-axis. 

The Type Bl solution is obtained by letting /(r, 9) — — \9 — f , a — |, b — — \, 
and n — 1. The profile functions of the gauge potentials are given by, 

1 11 
4>i = 2' A = sin6» - -cos -0(1 - cos0), 

1 1 

R x = 0, P 2 = cos0 + -sin -0(1 -cos 9), 

$1 = £sin^0, $ 2 = £cos^0. (38) 

The Type B2 solution is obtained by letting /(r, 0) = \9 + f , a = |, 6 = — |, and 
n = 1. The profile functions of the gauge potentials are given by, 

1 13 

■01 = — 7:, -Pi = sin + — cos —0(1 — cos 9), 

1 3 

R x = 0, P 2 = cos0 - -sin -0(1 - cos0), 

$1 = -Zsinh, $ 2 = -£ C os^0. (39) 

Both Type B solutions possess a magnetic charge of ^ at r = and 't Hooft's 

gauge potential, Ai = { 2° s s f ne I wn i cn is singular along the negative z-axis. 
The magnetic field B? of the above four solutions are given by 

gB? = (40) 

when the semi-infinite Dirac string along one-half of the z-axis is ignored. Thus 
the four solutions possess common magnetic field of a one-half monopole located at 
the origin, r = 0, that are spherically symmetric, radial in direction and pointing 
in the Higgs field direction in isospin space. 

Similar to the gauge potentials, the Higgs magnetic fields of the Type A solu- 
tions are singular along the negative z-axis, whereas the Higgs magnetic fields of 
the Type B solutions are singular along the positive 2-axis and they are given by 

„ \ sin(^0) fs , \ 

g Bf = 1 Type A solutions 41 

sm0 r 2 

gB? = l ^—4, Type B solutions. 42 

sm0 r 2 

The net 't Hooft's magnetic fields for the Type A and Type B solutions are however 
given by 

gB t = ^- 2nS(x 1 )S(x 2 )9(x 3 )Sf : Type A solutions, (43) 
9Bi = 2^2 + 2 ^(x 1 )6(x 2 )e(-x 3 )6l Type B solutions, (44) 

9 



when the Dirac string in the gauge potentials, Ai, is taken into consideration [T7] . 
These one-half monopole solutions therefore carry a positive one-half magnetic 
monopole at the origin and a semi- infinite Dirac string of flux, — , going into the 
origin. Hence the net magnetic charge of the configuration is zero. 



5 The Numerical Solutions 
5.1 The Numerical Calculations 

The numerical one-half monopole solutions can be constructed by making use of 



the exact one-half magnetic monopole solutions (36) - (39) as asymptotic solutions 



at large distances and by fixing the boundary conditions for the profile functions 



(16) along the z-axis and near r = 0. Hence, we are able to obtain four types 
of finite energy numerical one-half magnetic monopole solutions which we will 
name as the numerical Type Al, A2, Bl, and B2 one-half monopole solutions 
respectively. Since the function R 2 (r, 9) is singular along one side of the z-axis at 
infinity, we choose to perform our numerical analysis with the functions, 

=Va(r,0) sin 0, P 2 (r, 9) = R 2 (r, 9) sin0. (45) 

Near the origin, r = 0, we have the common trivial vacuum solution for the four 
solutions. The asymptotic solutions and boundary conditions at small distances 
that will give rise to finite energy solutions are 

fa = P 1 = R 1 = P 2 = Q, $ 1 = £ o cos0, $ 2 = -£ o sin0, (46) 

sin0$i(O,0) +cos0$ 2 (O,0) = 0, 
<9 r (cos0$i(r,0) -sin0$ 2 (r,0))| r=o = 0. (47) 

The boundary conditions imposed along the positive z-axis for the profile functions 



(16) of the Type A solutions are 



d e Mr,9)\ e=0 = 0, <£> 2 (r,0) = 0, d e ^{r, 9)\ e=0 = 0, 

i2i(r,0) = 0, P(r,0) = 0, P 2 (r,0) = 0. (48) 

Along the negative z-axis, the boundary conditions imposed are 

$ 1 (r,vr) = 0, d e $ 2 {r,9)\ e ^ = 0, d e ^i{r, 0)| e=7r = 0, 

R 1 (r,ir) = 0, P 1 (r,*) = 0, d e P 2 (r, 9)\ e=7r = 0. (49) 

The boundary conditions imposed along the positive z-axis for the profile functions 



(16) of the Type B solutions are 



$i(r,0)=0, «%$ 2 (r,0)| e=o = O, d^i(r, 9)\ 0=o = 0, 

R 1 (r,0) = 0, Px(r,0)=0, d e P 2 (r, 9)\ e=0 = 0, (50) 
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and along the negative z-axis, the boundary conditions imposed are 



00$i(r,0)| e=7r = O, $ 2 (r,7r)=0, ^(r, 0)| e=7r = 0, 

R 1 (r,n) = 0, P 1 (r,n) = 0, P 2 (r,7r)=0. (51) 

The numerical Type AL, A2, P2, and P2 solutions connecting the asymptotic 



solutions (36) - (39) respectively at large distances to the trivial vacuum solution 



(46) at small distances and subjected to the boundary conditions (47) - (49) for 



the Type A solutions and the boundary conditions (47), (50) - (51) for the Type 



B solutions together with the gauge fixing condition [8] 

rdrRx - d e ipi = 0, (52) 



are solved using the Maple 12 and MatLab R2009a softwares pi 

The second order equations of motion ^ which are reduced to six partial 



differential equations with the ansatz (16) are then transformed into a system 



of nonlinear equations using the finite difference approximation. The numerical 
calculations are performed with the system of nonlinear equations been discretized 
on a non-equidistant grid of size 90 x 80 covering the integration regions < x < 1 
and < 9 < it. Here x = is the finite interval compactified coordinate. The 
partial derivative with respect to the radial coordinate is then replaced accordingly 
by d r ->■ (1 - x) 2 d s and ^ -> (1 - x) 4 - 2(1 - x) 3 § . We first used Maple 
to find the Jacobian sparsity pattern for the system of nonlinear equations. After 
that we provide this information to Matlab to run the numerical computation. The 
system of nonlinear equations are then solved numerically using the trust-region- 
reflective algorithm by providing the solver with good initial guess. The second 
order equations of motion Eq. ^ are solved when the 0-winding number n = 1, 
with nonzero expectation value £ = 1, gauge coupling constant g — 1, and when 
the Higgs self-coupling constant < A < 12. In our calculations, there are two 
kinds of error, the first inherited from the finite difference approximation of the 
functions, which is of the order of 10 -4 . The second comes from the linearization 
of the nonlinear equations for MATLAB to solve numerically; the remainder or 
approximation error is of the order of 10~ 6 . Hence the overall error estimate is 
10" 4 . 

The numerical solutions obtained for the four one-half monopole solutions, 
■01, Pi, Ri, Pi, $i, and $2, are all regular functions of r and 9. However the function 
R 2 = possesses a string singularity along the negative (positive) z-axis for the 
Type A (Type B) solutions. Our results show that the Type Bl and Type B2 
solutions are exactly 180° rotation of the z-axis about the origin, r = 0, of the Type 
Al and Type A2 solutions respectively. The 3D profile function graphs versus r 
and 9 in rectangular coordinate frame of grid size 90 x 80 are shown for the Type 
Al and Type A2 solutions, Pi(r,9) and P 2 (r,9) in Figure [TJ ^(r, #) and R 1 (r,9) 
in Figure [2j and $i(r, 9) and $i(r, 6) in Figure [3j 
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5.2 The Magnetic Dipole Moment and Magnetic Charge 



From Maxwell electromagnetic theory, the 't Hooft's gauge potential of Eq. (28) 
at large r tends to 

Ai = (cos a + cos K)d i <l)\ r -+ 00 = f 1 _ ( ^(cos 9 ± 1) + F ° ® ) , (53) 

rsm# \2 r / 

F G (0) = r{A(Px -sm0) -s(P 2 - cosfl) - ^(cos0 ± I)}]™, (54) 

for the Type A and Type B solutions respectively From the graphs of Fg{9) versus 
angle 9, Figure|4](a), we find that Fg{9) = fi m sin 2 9, where \i m is the dimensionless 
magnetic dipole moment of the one-half monopole and is non vanishing for all 
values of A. The value of \i m is read from the graphs of F G {9) versus angle 9 at 
9 = |. The graphs of the magnetic dipole moment fi m i and fi m 2 of the Type 
A solutions versus A 1//2 when g — £ = 1 are shown in Figure [Z] (c) for values of 
< A < 12 and the numerical values are given in Table [TJ In the limit when A = 0, 
the dimensionless magnetic dipole moment for the Type 1 solution is [i m \ = ±1.32 
and for the Type 2 solution is fj, m2 = ±1.04. The plots of [i m \ and fi m2 versus A 1//2 
are non increasing graphs that concave upwards for both the Type A solutions. 
Also noted is that, // m i > // m 2, for the same value of A. Similar findings are 
also recorded in Ref. [8] and [18] for the MAP solutions where the dimensionless 
magnetic dipole moment decreases from A = (// m = 4.72,4.71 respectively) to 
A = 1 (fi m = 3.25,3.40 respectively) when the electic charge parameter rj = 0. 
At spatial infinity, these four one-half monopole solutions are exact and given 



by Eq. (36) to (39). From these exact solutions, we are able to conclude that 



the net magnetic charge of the four one-half monopole solutions is zero. This is 
because at spatial infinity, the Dirac string carries a magnetic flux of y going 
into the center of the sphere at infinity, r = 0, and the one-half monopole at the 
origin carries a magnetic flux of y coming out from the origin, r = 0. In order to 
check that our numerical calculations are consistent with the exact solutions, we 



calculate numerically and plot the graphs of the net magnetic charge M (|8j), the 
topological magnetic charge Mh (|9|), and the magnetic charge carried by the gauge 



field Mq (10) versus the compactified for the Type 1 and 2 solutions 

~% = 1. See Figure El (d). We are only able to obtain the magnetic 



when A 

charge, Mh, of the one-half monopole from the numerical data. The numerical 
data cannot account for the magnetic charge carried by the Dirac string. We also 
notice from Figure (d) that M — > as r — > and M — | when r > 4. 



From Eq. (|8j), and using the definition for the Abelian magnetic field (31 ) given 
by Bogomol'nyi j3] and Faddeev [IB], we can write 

M = — [ d% d 3 x= [ M d9 dr, where M = -r 2 sin 9{d i B i }, (55) 
47T J J 2 

is the magnetic charge density. The 3D surface and contour plots of A4 of (a) 
Type Al, (b) Type A2 solutions and (c) 't Hooft-Polyakov monopole the along 
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the x-z plane at y — when A = £ = 1 are shown in Figure [5] We note that the 
magnetic charge density is concentrated along the negative (positive) z-axis for 
the Type A (Type B) solutions. We also note that the magnetic charge density is 
solely positive for the four one-half monopole solutions. 



Type A I 



Type A2 
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Figure 1: The 3D profile function graphs Pi(r,6) and P2(r,8) versus r and 9 in 
rectangular coordinate frame of grid size 90 x 80 of the (a) Al and (b) A2 one-half 
monopole solutions when A = £ = 1. 



5.3 The Higgs Field and the Magnetic Field 



From Eq. (28), the 't Hooft's magnetic field lines contour plots of the Type Al 
and Type A2 one-half monopole solutions along the x-z plane at y = are shown 
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0.60 


0.61 
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Table 1: The magnetic dipole moments \i m xi I^m2 } and total energies Ex, E 2} of the 
Type 1 and 2 solutions respectively for < A < 12. 

in Figure [6] when A = £ = 1. In the finite r region near the origin, r = 0, 
the magnetic field lines of the numerical Type A (Type B) one-half monopole 
are string-like and concentrated along a finite length of the negative (positive) 
z-axis before they become hedgehog-like at large r. For the Type 1 solutions, the 
magnetic field line string stretches from z ~ ±0.14 to z ~ =f3.82 and for the Type 
2 solutions, the magnetic field line string stretches from z £s ±0.14 to z ~ =f3.25. 

The direction of the magnetic field lines of the one-half monopole solutions 
are shown by the vector field plots of the 't Hooft's magnetic field unit vector Bi 



* 



(29) along the x-z plane at y = when A = £ = 1. See Figure |6j The asterisk 
indicates the location of the one-half monopole at r = 0. The presence of a 
one-half monopole sitting at the origin is once again indicated in Figure [6] by the 
location of the Higgs field's node marked asterisk * on the Higgs field's vector, $ a , 
plots along the x-z plane at y = 0. 

The 3D surface and contour plots of the modulus of the Higgs field functions 
|$| of the (a) Type Al and (b) Type A2 one-half monopole solutions and (c) the 't 
Hooft-Polyakov monopole solution along the x-z plane at y — when A = £ = 1, 
show that there is a point zero of the Higgs field modulus at r = as shown in 
Figure [7| Hence the one-half monopoles and the one-monopole are located at r = 
where the Higgs field vanishes. The shape of the 3D surface plots for the one-half 
monopole in Figure [7] (a) and (b) is that of a flatten cone with its half-rounded 
vertex touching r = 0. The cones are flatten along the negative z-axis for the 
Type A one-half monopole. This is in contrast to the 3D surface plot of the Higgs 
modulus for the 't Hooft-Polyakov monopole which is a sharp vertex circular cone 
with its vertex at r = 0, where there is only one zero of |$| at r = 0. In the MAP 
solutions, there is also only one zero of |$| at the location of the MAP dipole when 
the </>- winding number, n — 1. However when n = 2, there is a double zero of |$| 
at the location of the MAP dipole along the z-axis in the direction perpendicular 
to the z-axis [IB]. In the case of the one-half monopole there is a double zero of 
|$| at r = along one side of the z-axis only, the negative (positive) z-axis for the 
Type A (Type B) solutions. 

The modulus of the Higgs field at large r tends to 

3^,9)1^ = ^-^1). (56) 
From the graphs of the Higgs function Fh{Q) versus the angle 9, Figured (b), we 
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find that Fh{0) is a non vanishing constant C\ only when A = for both the Type 
1 and Type 2 solutions. The constant c\ = —0.43 for the Type 1 solutions and 
c\ = —0.63 for the Type 2 solutions. For non vanishing values of A, the constant 
C\ vanishes for both Type 1 and Type 2 solutions. 

5.4 The Energy Density and Total Energy 



The total dimensionless energy (12) when £=g=l, can be written as 



E = J £d6 dr, £ = ^r 2 sin9{B?B? + D& a D& a + ^($ a $ a - £ 2 ) 2 }, (57) 

where £ is the energy density. The 3D surface and contour line plots of the energy 
density £ of the (a) Type Al, (b) Type A2 one-half monopole solutions and (c) 
the 't Hooft-Polyakov monopole solution along the x-z plane at y = when A = 1 
are shown in Figure |8j The Type 1 solutions has a slightly wider spread of the 
energy distribution along one-half of the z-axis and a slightly lower energy density 
peak with the same value of A = 1 compared to the Type 2 solutions. For the 
Type Al solution the energy density peak is located along the z-axis at —2.941 
with peak value of 1.087 while the energy density peak of the Type A2 solution 
is located along the z-axis at z = —2.235 with peak value of 1.149. The energy 
density concentration of the Type A (Type B) solutions are along the negative 
(positive) z-axis centered around its peak value. The 't Hooft-Polyakov monopole 
solution has two peak values of 0.3629 at z=0 and x = ±0.9412. Hence the shape 
of the one-half monopole is that of a rugby ball whereas the 't Hooft-Polyakov 
monopole is shaped like a torus. 



The total dimensionless energy (57) of the Type 1 one-half monopole solutions 
is Ei = 0.509 and that of the Type 2 solutions is E 2 = 0.525 when A = 0. Their 
total energies seem to be higher than the BPS total energy of a one-half monopole 
which is |. The total energies of the one-half monopole solutions are plotted 
versus \ 1 ' 2 and ln(l + A) for values of Higgs field strength < A < 12, Figure [9] 
(a) and (b). For a particular value of A, the Type 2 solutions possess higher total 
energy, compared to the Type 1 solutions, E\ > E 2 . The total energy E\ and E 2 
versus A 1 / 2 curves are non decreasing graphs that concave downwards, Figure [9] 
(a), whereas the total energy E\ and E 2 versus ln(l + A) graphs become linear 
when A is large, Figure [9] (b). The table of total energies Ei and E 2 for different 
values of A of the one-half monopole solutions are shown in Table [TJ 

5.5 Comparison between the Numerical Solutions 

The boundary condition at r — > 00 is fixed by the exact one-half monopole solu- 



tions (36) - (39) and the four exact solutions are related by gauge transformations. 
However near the origin as r — > and along the positive and negative z-axis (or at 
6 = and 7r respectively), the solutions are fixed only by similar boundary condi- 
tions. A direct method of checking the similarities between the Type A and Type 
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B solutions is to rotate the z-axis of the Type A solution about its center at r = 0, 
180° so the the positive z-axis becomes the negative z-axis. The transformation 
performed on the Type A profile functions, P\, ipi, and $1, to make it similar to 
the respective Type B profile functions is 

F A (r, 9) -> F A (r, tt - 9) = F A (r, 9) w F B (r, 9), (58) 

and the transformation performed on the Type A profile functions, P2, R%, and 
$2, to make it similar to the respective Type B profile functions is 

F A (r, 9) -> -F A (r, ir-9) = F A (r, 9) « F B (r, 9), (59) 

where F A (r,9) and F B (r,9) represent the Type A and Type B profile functions 
respectively. The transformed Type Al(2) profile functions, F A (r,9) are then 
subtracted from the respective Type 51(2) profile functions, F B (r,9). In our 
numerical calculations, we found that the 3D surface plots of the differences, 
F B (r,9) — F A (r,9), versus r and 9 is very small. They are of the order of 1CT 13 
for the Type 1 solutions and of the order of 10 -12 for the Type 2 solutions. Hence 
the Type B solutions are exactly 180° rotations of the z-axis about r = of the 
respective Type A solutions. 

In order to compare the four types of solutions with each other to check if the 
Type 1 and Type 2 solutions are genuinely different and also to confirm that the 
Type A and Type B solutions are gauge equivalent, we used the gauge invariant 
form of the gauge potential and the Higgs field, 

A 2 = r 2 sin 2 9A a i A a i1 $ 2 = $ a $ a . (60) 

From the four types of solutions, there are six combinations of the differences 
between the solutions of the quantities A 2 and $ 2 . We then plot twelve 3D surface 
graphs of the differences of A 2 and $ 2 versus r and 9 for all the six combinations. 
We find that for the graphs of (A A1 - A Bl ), (& A1 - $|i), (A A2 - A 2 B2 ), and 
($^ 2 — $52), versus r and 9, the order of magnitudes are very small, 10 -12 , 10~ 13 , 
10~ 15 , and 10~ 15 respectively. Therefore the Type A and Type B solutions are 
gauge equivalent solutions and they differ by only a rotation of 180° of the z-axis 
about its center in 3D space. The transformed functions A 2 and $ 2 obey the 



transformation (58). 



However for the 3D graphs of (A 2 A2 - A 2 B1 ), - & B1 ), (A 2 A1 - A 2 B2 ), {$ 2 A1 - 
&B2), (-^42 - A 2 A1 ), ($| 2 - & M ), (A 2 B2 - A 2 B1 ), and ($| 2 - versus r and 9, 
there are significant differences between the Type 1 and Type 2 solutions along the 
z-axis in the region where the energy of the one-half monopole is concentrated. As 
expected the four sets of graphs for the differences in A 2 and $ 2 are similar. We 
therefore show in Figure [9] (c) and (d) only one of the four sets, the 3D surface plot 
of {A A2 — A 2 Al ) and (& A2 — $ Al ) when A = 1. The maximum value of (A 2 A2 ~ A 2 A1 ) 
is 0.1825 and {§ 2 A2 — $ai) is 0.2750. Both maxima occur at the same point along 
one side of the z-axis at r = 3.38 and 9 = tt. Hence the Type 1 and Type 2 
solutions are not gauge equivalent one-half monopole. They are different and the 
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differences in the solutions occur in the region of where the energy of the solutions 
is concentrated. 

We also note that the gauge transformations (19) connecting the Type Al (36) 
to Bl (|38~| solutions and the gauge transformation connecting the Type A2 (37) 



to B2 (39) solutions are. 



7T _ 7T 

(Al-KBl) — — 7^' J(A2^B2) — 2' 



(61) 



The gauge transformations (19) connecting the Type Al (36) to B2 (39) solutions 
and that connecting the Type Bl (38) to A2 (37) solutions are 



/(A1-KB2) — f{Bl^A2) — 



71 



(62) 



The gauge transformations (19) connecting the Type Al (36) to A2 (37) solutions 
and that connecting the Type Bl (38) to B2 (39) solutions are 



f{Al-*A2) — f{Bl~>B2) = + 71. 



(63) 



The gauge transformations (61) are trivial as / is a constant whereas the gauge 
transformations (62) and (63) are non-trivial and depend on 9. 



6 Comments 

We have found four types of one-half monopole solutions, Type Al, Type A2, Type 
Bl, and Type B2 solutions. All these four solutions possess zero net magnetic 
charge. By comparing all the six profile functions of the Type A and Type B 
solutions, we are able to conclude that the Type B solutions are exact 180° rotation 
of the 2-axis about r = of the respective Type A solutions in 3D space. Again 
by comparing the gauge invariant form of the gauge potentials, A 2 , and the Higgs 
fields, $ 2 , we are able to conclude that the Type A and Type B solutions are 
gauge equivalent. However the Type 1 and Type 2 solutions are significantly 
distinct in the region where the energy of the one-half monopole is concentrated. 
The difference in the Type 1 and Type 2 solutions are also reflected in the graphs of 
magnetic dipole moments versus A 1//2 , Figure [4] (c); and the graphs of total energy 
versus A 1//2 and ln(l + A), Figure [9] (a) and (b). All the one-half monopole solutions 
reported here satisfy the gauge condition, rd r Ri—dgipi = [8]. However this gauge 
condition still do not ensure the uniqueness of solutions for these one-half monopole 
configurations. We have in fact found two distinct one-half monopole, the Type 
1 and Type 2 one-half monopole, which are connected by the non trivial gauge 



transformations (62) and (63) 



For further work on this subject, we would like to introduce electric charge into 
the Type 1 and Type 2 one-half monopole solutions using the standard procedure 
of Julia and Zee [19J by letting the time component of the gauge potential at large 
distances be nonzero and introducing a constant electric charge parameter into 
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the exact asymptotic solutions at large r. The numerical finite energy rotating 
one-half dyon solutions will be presented in a separate work soon. 

Another direction of our further work is to look for a one-monopole and a 
one-half monopole of opposite magnetic charge in the SU(2) YMH theory. We 
have already found these type of solutions and will be finalizing our findings in a 
separate work soon. 
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Figure 2: The 3D profile function graphs ipi(r,9) and Ri(r,8) versus r and 9 in 
rectangular coordinate frame of grid size 90 x 80 of the (a) Al and (b) A2 one-half 
monopole solutions when A = £ = 1. 
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Figure 3: The 3D profile function graphs $i(r, 9) and $ 2 (r, #) versus r and 9 in 
rectangular coordinate frame of grid size 90 x 80 of the (a) Al and (b) A2 one-half 
monopole solutions when A = £ = 1. 
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Figure 4: Plots of (a) F G {6) and (b) F H {6) versus 5 when, A = 0, 1. (c) Plots of 
/i m i and /i m2 versus A 1//2 . (d) Plots of M, M G and M H versus x when A = 1. 
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Figure 5: 3D surface and contour plots of the magnetic charge density, M, of the 
(a) Type Al, (b) Type A2 and (c) the 't Hooft-Polyakov monopole solutions along 
the x-z plane at y — when A = £ = 1. 
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Figure 6: Contour plots of the 't Hooft magnetic field lines, unit vector field plots 

of the 't Hooft magnetic field and Higgs field vector plots of the (a) Type v41 and 

(b) Type solutions along the x-z plane at y = when A = £ = 1. The * 

indicates the location of the one-half monopole. 
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Figure 7: The 3D surface and contour plots of the Higgs field modulus |$| of the 
(a) Type Al, (b) Type A2 and (c) the 't Hooft-Polyakov monopole solutions along 
the x-z plane at y — when A = £ = 1. 
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Figure 8: 3D surface and contour plots of the energy density £ of the (a) Type 
Al, (b) Type A2 and (c) the 't Hooft-Polyakov monopole solutions along the x-z 
plane at y — when A = £ = 1. 
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